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We investigate the two-dimensional cooperon-fermion model in the correlated regime with a new 
continuous-time diagrammatic determinant quantum Monte Carlo (DDQMC) algorithm. We esti- 
mate the transition temperature Tc, examine the effectively reduced band gap and cooperon mass, 
and find that delocalization of the cooperons enhances the diamagnetism. When applied to dia- 
magnetism of the pseudogap phase in high-Tc cuprates, we obtain results in a qualitative agreement 
with recent torque magnetization measurements. 
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I. INTRODUCTION 

The coopcron-fcrmion model (for a review, see Ref. 
[l|) is a basic model for superconductivity that has 
widely been adopted to explain the BCS-BEC crossover 
in ultra-cold fermionic atomic gasesi"— and high-Tc 
superconductivityjii^^— The resonantly paired fermions, 
or cooperons, in this model can either be locally bound 
pairs of small polarons due to extremely strong electron- 
phonon coupling^ or localized Cooper pairs due to 
strong local pairing as might be the case in high-Tc 
superconductors^^, or molecular bosons in ultra-cold 
atomsi^"— The potential existence of finite energy cooper- 
ons with a local attraction has also been put forward 
a few years ago in a simple semiconductor systemJ^ 
One recent work based on the four-leg ladder Hubbard 
modelii as well as earlier studies in the quasi-2D lad- 
der Hubbard modeliS reveals the important role of the 
cooperon excitations in the transition from insulating 
state to superconducting state. 

In cuprates, the interplay between the finite energy 
cooperon excitations around the antinode and the trun- 
cated Fermi surface around the node has recently been 
proposed as a possible mechanism for the superconduc- 
tivity in a two-gap scenario J^i^ The dominant under- 
lying mechanism for driving superconductivity in this 
model is the scattering two electrons involving a vir- 
tual cooperon c-|-.fc -f- c^.-fc '''°"°> c-i- k> + c^.^k'- Com- 
pared with the attractive Hubbard model, the cooperon- 
fermion model has much richer physics since it has the 
complete dynamical information of the interaction be- 
tween two fermions and the delocalization of cooperons 
with decreasing temperature. 

So far most work on this model is done at the mean 
field level, with either the T-matrix method or various 
other methods going beyond simple mean field theory by 
including more diagrams^ but there have been few unbi- 
ased calculations. A recent exact diagonalization study 
of this model has been limited by small sizes and spe- 



cial geometries)^ A direct quantum Monte Carlo simula- 
tion usually suffers from a sign problem;^ except for cer- 
tain models and algorithms, such as determinant quan- 
tum Monte Carlo simulations of the attractive Hubbard 
model. 

In this paper, we simulate the two-dimensional 
cooperon-fermion model using a continuous-time dia- 
grammatic determinant quantum Monte Carlo method 
(DDQMC). The Hamiltonian has the form: 
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t 



bi + h.c.) 
(1) 

where cj^(ccr) is the fermionic creation (annihilation) op- 
erator with spin a = {t, 4-} and (b) is the bosonic cre- 
ation (annihilation) operator of a cooperon. The interac- 
tion term U leads to the s-wave pairing of fermions me- 
diated by the originally localized cooperons with a band 
gap A and the delocalization of cooperons at low temper- 
atures. The bare dispersions are e^, = 2t/[2 — cos{kx) — 
cos(fcy)]-|-/x, and = 2ti,{2 — cos{kx) ~ cos{ky)) + A. It is 
a trivial generalization to include attractive interactions 
between the fermions. 

By integrating out the bosonic degrees of freedom, an 
effective action for the fermionic part can be obtained, 
which has a form similar to the attractive Hubbard model 
but with full dynamic properties: 

rP 

Jo Jo 

(2) 



where G^'°_^.^ is the bare bosonic Green's function. 
The attractive Hubbard model can be obtained by ap- 
proximating C/^G^'°_^^^^_^, = UhSi.hSt^ti with Uh ^ 
— [/^/Aoff (Acff is the renormalized band gap of the 
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FIG. 1: (Color online.) (a): elementary interaction events 
(vertices): two fermions are combined onto a cooperon (left) 
and vice versa (right), (b): a diagram contributing to the 
partition function Z. (c): a diagram sampled in the DDQMC 
algorithm: open ends of the fermionic lines imply a sum over 
all the ways of connecting the vertices by the fermionic lines, 
which is represented by the corresponding determinant. 



cooperons) including the dominant contribution from the 
rcnormalized fe = cooperon. 

We develop a continuous-time DDQMC algorithm, 
similar to the algorithms for an attractive Hubbard 
mode]32.i2i for the cooperon-fermion model, and obtain 
the numerically exact solution to the model at the filling 
value n = 0.12. The Kosterlitz-Thouless (KT) transition 
temperature is estimated from the finite-size scaling of 
the fermion pair correlation function and the cooperon 
Green's function. The renormalization of the cooperon 
band characterized by the effective gap Aes and the ef- 
fective mass TTT-eff are examined carefully. Applying these 
results to study the strong diamagnetism recently ob- 
served in cuprates, we find that the renormalization of 
the cooperon band will enhance the diamagnetism dra- 
matically at low temperatures, which qualitatively agrees 
with the experimental data<2^ 



II. THE ALGORITHM 



Continuous-time DDQMC algorithms have been ap- 
plied successfully in the past to obtain the critical tem- 
perature in the BCS-BEC crossove r^^'^^ The sign prob- 
lem in these algorithms is avoided by collecting all Feyn- 
man diagrams with the same distribution of vertices as a 
single configuration which turns out to have a positive- 
definitive weight. This approach can also be applied to 
the cooperon-fermion model eliminating the severe sign 
problem coming from the permutations of the fermionic 
lines. In the interaction picture, the partition function 



for the model ([T]) can be expressed as 



Z = Tr 



(3) 



where iJ^? ^ in the Hamiltonian of free fermions with the 
spin cr and is the Hamiltonian for free cooperons, 
the bilinear in Cj[. ^, Cfc^o- and 6^, 6^ terms respectively in 
Eq. IT]). To describe the Feynman diagrams, we define 
two different kinds of events, shown in Fig. [IJa), repre- 
senting the process of combining two fermions with op- 
posite spins into one cooperon and the reverse process. 
For a typical Feynman diagram like the one shown in 
Fig. [TJb), the vertices are connected by the bare single 
particle propagators 
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The superscript a/c on r^j- are for the events with anni- 
hilation and creation of a pair of fermions (in accompany 
of the creation and annihilation of a cooperon), respec- 
tively. By applying the Wick's theorem, the partition 
function can be rewritten as 
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where S'„ represents the configuration including all pos- 
sible ways of connecting a specific distribution of ver- 
tices with the propagator lines. The matrix com- 



ponents [j4^' ^] 



G 



and [Bg' ] 



— G^^'iL^o ^a_^a. The determinant of the matrix A comes 

i j ^ i j 

from the anti-commutation relation of fermions, while 
the permanent Perm(i?g' ) originates in the commuta- 
tion relation between cooperons. In contrast to determi- 
nants, which can be efhciently evaluated, the calculation 
of a permanent is an exponentially hard problem. Thus, 
we evaluate the permanent by individually sampling all 
permutations of the bosonic lines. A typical diagram 
Sn encountered in the Monte Carlo sampling is shown 
in Fig. [BJc). The open ends of fermionic lines indicate 
that all possible connection ways of the fermion lines are 
summed up in the determinant, while the connection be- 
tween the cooperon lines is fixed, indicating that the spe- 
cific connections are sampled individually. Summing the 
fermion lines into a determinant completely removes the 
fermionic sign problem. However, sampling the perma- 
nent gives us a small remaining sign problem, which is 
tractable since the distribution of values in Bg' is domi- 
nantly positive. 
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The weight of a configuration Sn is 

U2n 



dctAg^^^dctA^^^^ W Bf\ (7) 
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where is the spatial size of the system with periodic 
boundary conditions. Thermal averages are calculated 
by sampling all possible configurations Sn- 
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FIG. 3: (Color online.) The Feynman diagrams comprising 
the fermionic (a) and cooperon (b) Green's function within 
the random phase approximation (RPA) . 




FIG. 2: (Color online.) Two sets of complementary Monte 
Carlo updates: (a) adding or removing one pair of vertices, 
and (b) swapping the end points of cooperon lines. 

Our algorithm implements three different Monte Carlo 
updates: creating or deleting one pair of vertices which 
changes the order from n to n ± 1, and swapping the 
connection of cooperon lines as shown in Fig. [51 In order 
to improve the efficiency of sampling we pick a pair of 
times with a probability proportional to the cooperon's 
bare Green's function G^' r- Adding a pair of vertices to 
go from configuration Sn to Sn+i is accepted with an 
acceptance ratio of 



-Rarid = min [IIJ 
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r<b.o 



(8) 



and a corresponding equation for the removal. 

For the self-complementary process of swapping the 
connection of cooperon Greens functions the acceptance 
ratio is 



i?flip = min 1, 



(9) 



The fermionic and cooperon Green's functions can be 
measured as 



is an (n-|- 1) X (n-|- 1) matrix extending 
A by an extra column and row corresponding to the open 



The matrix Ag 



vertices c 



and 



"(_R' r') The notation (..)mc denotes 

the Monte Carlo average. The particle-particle correla- 
tion function of fermions {TTCfi^r,iCji,T,tc\^i ^.i -^c^/ ^.i ^) is 
measured as 
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(12) 



MC 



with the double occupancy characterizing the local pair- 
ing strength (n^n;) =Xfl'L(o,o),T=o- • ^^'^ vertex cor- 
relation contribution to the particle-particle correlation 
of fermions, which indicates the formation of coherent 
cooper pairs, is defined as 



od — T t t 



(13) 



The cooperon Green's function and x°'^ in lo^ig wave- 
length and static limit for our finite system with peri- 
odic boundary conditions can be defined as G^^q — 

lo dr j:rG'r,r, and xlU^=, = Jo dr^r X^r with r be- 
ing the distance confined by system size. 

In addition to the DDQMC results, we also show re- 
sults of random phase approximation (RPA) calculations 
which only take into account the contribution from a 
truncated set of Feynman diagrams, whose diagrammatic 
representations are shown in Fig. [3l The cooperon and 
fermionic Green's functions have the general form 



h.uj 



(14) 



with the self energies S^^^ estimated by only including 
ladder diagrams. RPA works well at weak coupling and 
high temperature region, and is useful as a test for our 
numerical solution in that limit. 



G^R-R'^r-r'.a = ( ,,/det A^^^ 
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III. MAIN RESULTS 

For our simulations we choose 1/ = 1 as the unit of 
energy and set — 1, = 0.5. The bare cooperon 
band lies above the bottom of the fermionic band with 
the offset A = 0.75 and the total charge density is fixed 
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FIG. 4: (Color online.) (a): the chemical potential versus 
temperature T /U . (b): the particle density 71^ of a single spin 
component and the cooperon density versus temperature 
T /U . DDQMC results in thermodynamic limit (blue points) 
are extrapolated to the infinite system size. 



at n| + + 2nb — 0.12 ± 0.002. In the parameter 
regime we are interested in, the chemical potential /i 
is around 0.22 — 0.45, resulting in the effective renor- 
malizcd gap Aoff < 0.3 and the corresponding effective 
attractive Hubbard interaction \Uh\ ^ t/^/Aeff > 3 at 
/? > 3. For the finite-size scaling analysis, we use the set 
of linear system sizes L = 11, 15, 21, 25. The expectation 
values and the error bars are obtained from 96 indepen- 
dent sampling processes with different random number 
seeds. Each measurement is made after 2-3 autocorre- 
lation times, i.e. around one measurement per lOOL'^ 
Monte Carlo steps at high temperatures (/3 ^ 3), and 
3000-5000 steps for low temperatures (/3 - 17). 

Fig. mja) shows the temperature dependence of the 
chemical potential yu. The solid curve is the RPA re- 
sult with a mean-field critical temperature T^^^ ^ 0.09 
characterized by the closure of the effective gap Aoff. 
Blue diamonds are the chemical potential at thermody- 
namic limit determined by DDQMC. RPA generally over- 
estimates the interplay between fermions and cooperons 
at low temperatures due to the logarithmic divergence 
leading to a larger value of the cooperon's self-energy 
from the simple particle-particle bubble diagram. As a 
consequence, the RPA chemical potential is lower than 
the exact value in this regime. However, in the high- 
temperature limit, the RPA calculation provides a con- 
sistency check for DDQMC and there we observe a per- 
fect agreement. Fig. IDJb) shows the particle density 
and n(, at thermodynamic limit. The cease of the sup- 
pression of Ub as the temperature is decreased, i.e. the 
flattening out of nb v.s. T at low r(< 0.1), indicates an 
approach to the KT transition. 

As seen from Fig. El the double occupancy = 
{ni^-fUi^l) ^ reveals a strong on-site pairing 

mediated by cooperons. The behavior of is determined 
by two aspects: (i) the competition between the potential 
energy gain from pairing and the corresponding kinetic 



FIG. 5: (Color online.) The system size dependence of the 
onsite double occupancy Ud for various values of /3. 



energy loss, and (ii) the balance between the particle 
number of fermions and cooperons. However, the latter 
is not expected to play an important role at low tem- 
peratures due to the plateau in nb{T). The monotonous 
increase of Ud with lowering the temperature is analo- 
gous to recent DMFT results for the attractive Hubbard 
modelj^ We also note that a different low-T behavior of 
rid for the attractive Hubbard model has been reported 
in early QMC studies<^ 

The 2D superconducting state is characterized by al- 
gebraically decaying off-diagonal order— in x"''- Figure 
inija) shows Xijl uj=Q for different system sizes, where 
-Rx.max = — l)/2 is the maximum distance in the x di- 
rection on the lattice with periodic boundary conditions. 
We also plot the size dependence of the cooperon Green's 
function G^^ in Fig. Hl^b), which behaves very 

similarly to x'r uj=o- Both grow substantially with 
decreasing the temperature, indicating that the Cooper 
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FIG. 6: (Color online.) (a): system-size dependence of off- 
diagonal order XhI „,ax '^^ distance ii^^.max. (b): the 
cooperon Green's function GJ;^ ^^^^ ^^^q. i?i^,max = {L — l)/2 
is the maximal distance along the x direction in our system 
with periodic boundary conditions. 



5 



pairs and coopcrons become coherent. 
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FIG. 7: (Color online.) Finite-size scaling of Xfe=o u=o ^^'^ 
Gfe_Q jj^g according to Eq. (|15p : the error bars are smaller 
than the symbol size. The uncertainties of ^1 = 0.55 ±0.1 
and Tc = 0.03 ± 0.01 are estimated from a breakdown of the 
shown data collapse. 

We next perform a finite-size scaling analysis for the 
pair correlator XkLo uj=o ^'^'^ cooperon Green's func- 



For Tc > T > 0, one expects Xr. 



od 



ui=0 



to 



exhibit a power-law decay with an exponent ti{T), such 
that r]{Tc) = 1/4 and r]{T = 0) = indicating the emer- 
gence of the true log-range order at T = 0. Above Tc, 
Xr'L=o shows an exponential decay. The pair correlator 
at T > Tc should obey the scaling formulap2^ 

L2-"(^=)/(i/e^), for L > LT ^ T+ (15) 



X 



od 



with ^-^ ~ ^Al^jT~T^ ^ Since the critical behavior in both 
subsystems of cooperons and fermions is a manifestation 
of one and the same supcrfluid transition, G^^q is 



also supposed to exhibit the scaling given by Eq. ([151) . 
The parameters A and Tc are chosen so that measure- 
ments of Xfelo w=o ^'^'i G^^o for different system 
sizes collapse in the vicinity of the phase transition as 
shown in Figs. [71 resulting in A = 0.55 ± 0.1 and 
Tc = 0.03 ± 0.01. The uncertainties in A and Tc are 
estimated from observing a noticeable distortion of the 
data from a single smooth curve as the parameters are 
varied beyond the claimed error bars. 

The most important properties of cooperons are their 
effective band gap Aoff and effective mass TOcff renor- 
malized by interactions mediated by fermions. As we 
shall discuss in more details in the next section, these 
parameters will allow us to obtain an estimate of the 
diamagnetic susceptibility, which is expected to rise dra- 
matically due to the quasi-condensation on approach to 
Tc- Figure [UJa) shows the dependence of Acff, which is 
obtained from the cooperon Green's function according 
to Acff = — [G^^o tj=o]~"^' t'^^ linear system size. The 
value Acff.L_j.oo in the thermodynamic limit, obtained 
from the extrapolation in the system size, is shown in 
Fig. [5]^b). At T > 0.1, the gap Acff,L_i.oo is quite close to 



the RPA estimate. Due to the logarithmic divergence of 
the bare particle-particle bubble, RPA leads to a substan- 
tially higher mean-field critical temperature Tc ^^ ~ 0.09. 
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FIG. 8: (a) (Color online.) Finite size extrapolation of the 
renormalized cooperon gap Acfi. (b) Comparison of the ex- 
trapolated Acfi at L ^ cxj with the RPA results. Tc = 0.03 
obtained in Fig. [71 is shown by the arrow. The error bars are 
smaller than symbol sizes. 



The effective cooperon mass is calculated as 



-1 

off 



Efth 
r r,: 



[Gfc=0,tJ=o]^ 



(16) 



where the factor of in the sum shows the importance 
of the long-range behaviour of the Green's function in 
determining mcff . Since the finite size of the system along 
with the periodic boundary conditions will enhance ^ 
at large r, the straightforward evaluation of mcff using 
Eq. ([rS]) is inadequate. To get rid of the finite-size effects, 
we fit the measured ^ according to 
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(17) 



with = ae^''!''! for large |r| (both a and b depend 
on r) . The function Q'^ ^ obtained thereby is then used 
instead of Gf, in Eq. p^ . The result obtained using 
the data for i = 11, 15 is shown in Fig. MX>). The RPA 
calculation gives similar values at T > 0.15. At lower 
temperatures, the effective mass continuously decreases 
and tends to a finite value. Fits to Eq. ([TSI allow us to es- 
timate the renormalized dispersion of the cooperon band 
eff shown in Fig. [SI using the Green's functions for 
the case L = 11. For comparison, the RPA curve is 
also shown in Fig. [HI which agrees with the Monte Carlo 
results at high temperatures. 
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FIG. 9: (Color online.) (a): The suppression of the renor- 
malized cooperon mass with decreasing the temperature. Red 
curve corresponds to the RPA result. The blue (black) dots 
are for lattice size L — 11, 15. The same fitting process 
fails for larger lattice sizes L — 21, 25 due to the too small 
measured close to the lattice boundary and very large 
sampling error bars, (b): Comparison of the renormalized 
cooperon's dispersion £k,eff from RPA calculation and our 
simulation with L = 11. The fitting process (Eq. <\17^ ) is ap- 
plied to obtain ^k.eff in our simulation. The error bars are 
smaller than the symbol. 



IV. APPLICATION TO THE STRONG 
DIAMAGNETISM IN THE CUPRATES 

Lots of anomalous properties of the pseudogap phase 
have been reported since the early stage of high- Tc stud- 
ies, for instance the existence of the partial gap, the linear 
resistivity, and the proportionahty of the charge carrier 
density to doping concentrationi^2i2i Strong supercon- 
ducting fluctuations have been observed in a large tem- 
perature region in recent Nernst and torque magnetome- 
try measurements 1^1^ In contrast to conventional BCS 
superconductors, where the Gaussian (amplitude) fluctu- 
ations arc dominant and the pairing length is quite long, 
in the cuprates the fluctuations of the phase rigidity arc 
predominant while the cooper pairs are strongly bound 
with the energy scale around the spin-spin superexchange 
J ^ 100 meV and they localized in a small spatial area 
with ^ ~ 3 — 4 lattice constant (a = 3.8 A). 

In momentum space the pseudogap phase is highly 
anisotropic. More and more evidence shows that 
the states at the antinode and node are intrinsically 
different^. In the pseudogap phase, the single parti- 
cle gap is partially opened only around the antinode, 
leaving cither "arcs" or hole-like Fermi pockets around 
the nodesi^"— Recent angular-resolved photoemission 
spectroscopy (ARPES) experiments have observed the 
existence of such pocket and unmask the particle-hole 
symmetry of the spectrum around the antinode, and 
the particle-hole asymmetry around the node.— Evi- 
dence from scanning tunneling microscopy (STM) re- 
veals that the low-energy states around the nodes are 
homogenous and of long range correlation length, and 



quantum interference STM observed well-defined Fermi 
surface only inside AF reduce Brillouin zone^. Mean- 
while the high-energy states around the antinodes are in- 
homogeneous with short correlation length around four 
lattice constants. Some evidence of particle-hole sym- 
metry in ARPES^ around the antinodes is consistent 
with the functional renormalization group calculations^ 
showing that the strong umklapp scattering enhances the 
cooperon channel. 

Some possible scenarios leading to superconductivity 
by the finite energy cooperon excitations and fermion 
sea have been proposed from the Hubbard mode on a 
ladder-iii^ and semiconductors^^ which may shed some 
light on the case for cupratesi^ So far the possible in- 
terplay between the states residing on the node and the 
antinode is still an open question. The superconducting 
gap on the Fermi surface around the node in SC state 
may be driven by this effect. 




FIG. 10: (Color online) Anisotropy in momentum space in 
the pseudogap phase: localized tightly bound cooperons are 
located at the antinodes and a hole-like Fermi sea resides on 
the nodes. 

It is reasonable to treat the states around the node as 
a free Fermi gas and the component around the antinode 
as tightly bound cooperons with a pairing gap of the or- 
der of J ~ 100 meV as sketched in Fig[TOl Our model ^ 
can thus be a simplified picture for the cuprate pseudo- 
gap phase by ignoring the multi-patch structure and the 
d-wave phase of the momentum-space localized cooper- 
ons. There has been a lot of previous work attempting 
to use this model for the phenomenology of the pseudo- 
gap phase (for a review see Ref. HI). However most of 
these are based on mean field calculations, but an ex- 
act solution is still missing. Our numerical simulations 
can give useful quantitative insight. Here we will focus 
on the strong diamagnetism observed recently. From the 
recent phenomenological YRZ theoryji^ it is reasonable 
to further propose that the fermionic particles around 
the node are itinerant holes, and the cooperons around 
the antinode are tightly bound hole pairs. To describe 
the undcrdoped phase we assume the total charge carrier 
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density is around 0.12, which is obtained by a value of 
/i ~ lOOmeV (with U ^ 250meV), which is in a reason- 
able regime based on early analysis of the YRZ model^. 
The bare hopping tf ~ 250 meV is chosen comparable to 
the nearest neighbor hopping integral in cuprates. 

In conventional BCS theorj4i the contributions of the 
fcrmionic pairing fluctuation to diamagnetism are sub- 
stantial only in a very narrow temperature region above 
Tc- In our model, however, the renormalized cooperon 
band contributes dominantly to the diamagnetism in a 
wide temperature range above T^, in pseudogap phase. 
Since it is computationally extremely demanding to cal- 
culate the numeric value of the second order coefficient in 
q of the current-current correlation function ICq^ui=0: we 
approximate the diamagnetization of cooperons as that 
of free bosons with renormalized gap Acff in the limit 
L ^ oo and the effective mass nics at L = 11 (see Figl9]). 
In CGS units [l/(47r)] it has the formii^ 



X 



{2ef 



'fc=0 



^TOoffd 6 



(18) 



where d=6A is the interlayer distance for cuprates. The 
value (2e)^/(c^2me(i) corresponds to M—7.5 A/m at 
H = 1T. 

The experimentally observed singular behavior 
M{T,H) ~ -iJi/<5(T) ^iti^ ^ ^ at small H and 
T — !> T+ might be related to the mesoscopic Meissner 
effect or the fragile Landon rigidityi^ Thermally excited 
vortices with exponentially increased inter-vortex length, 
however, are not sufficient to explain the experimental 
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FIG. 11: (Color online.) The temperature dependence of the 
magnetization at B = 1 T calculated using Eq. (|18|) .The x- 
axis is normalized by Tc ~ 0.03 (~ 80 K with U = 250 meV). 
Red dots are from the torque magnetization experiments with 
Tc = 50K.— Note that there is no parametric fitting to the 
experimental data. While strong damping of cooperons at 
high temperature suppresses the magnetization at T 2> Tc, 
the exponential increase close to Tc will not be affected. 



observation. In Fig. [TT] we show the diamagnetism 
aX H = IT with M = xH and x calculated from 
Eq. (fTS]) . Strong diamagnetism prevails in a very wide 
temperature region above Tc] this is in agreement with 
the experimental data (with Tc = bQK). Note that even 
though the experimental data shows a much narrower 
temperature region, it is still orders of magnitude wider 
than that predicted by the conventional BCS theory. 
In the cooperon-fermion model, the strong damping 
of the cooperon at high temperature will suppress the 
diamagnetism greatly but will leave the exponential 
increase of diamagnetism at T close to Tc unchanged. 
We also note that an alternative explanation of the 
strong diamagnetism based on a vortex liquid picture 
has been proposed by Oganesyan et al.^ 



SUMMARY 



We have developed a continuous-time diagrammatic 
determinant quantum Monte Carlo algorithm for the 
cooperon-fermion model. Our results for the fcrmionic 
part of the model show similar behavior to its twin model, 
the attractive Hubbard model, which is often used to de- 
scribe the BCS-BEC crossover in the systems of ultra- 
cold atoms, where the cooperon-fermion model is the rel- 
evant model on the BEC side. 

Besides the critical temperature we have calculated the 
renormalized band gap and mass of the cooperons. The 
decrease of the mass and the suppression of the renormal- 
ized gap have important effects on the thermodynamic 
properties of the cooperons. Applied to cuprate super- 
conductors, the interplay between the cooperons at the 
antinode and the fermions at the node is expected to de- 
localize the cooperons and finally lead to a substantial 
enhancement of the diamagnetism in a wide temperature 
range. That could explain the strong diamagnetic signal 
observed recently in the underdoped state. 

The numerical method developed here can be used to 
study the BCS-BEC crossover on lattices in the frame- 
work of the cooperon-fermion model, which gives direct 
access to the paring physics via the cooperon part. The 
universal results in terms of to the s-wave effective cou- 
pling between the fermions can in principle be obtained 
in the low-density limit, as was done, e.g., in Refi^. 
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ory Center of the University of Maryland. We used the 
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